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Abstract 

We show that the observed evolution of the transmission phase through multi-electron quantum 
dots with more than ~ 10 electrons, which shows a universal (i.e., independent of N) as yet 
unexplained behavior, is consistent with an electrostatic model, where electron-electron interaction 
is described by a mean-field approach. Moreover, we perform exact calculations for an open 
ID quantum dot and show that carrier correlations may give rise to a non-universal (i.e., N- 
dependent) behavior of the transmission phase, ensuing from Fano resonances, which is consistent 
with experiments with a few (N < 10) carriers. Our results suggest that in the universal regime 
the coherent transmission takes place through a single level while in the few-particle regime the 
correlated scattering state is determined by the number of bound particles. 

PACS numbers: 73.63.Kv, 03.65.Nk, 72.10.-d 73.23.Hk 
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I. INTRODUCTION 



Among the experiments that exploit the coherent dynamics of carriers, the one performed 
in 1997 by Schuster et al.— , in which the transmission phase of an electron scattered through 
a quantum dot (QD) was measured, constitute an ideal test on the validity of different 
theoretical models for the inclusion of electron-electron interaction. In fact, the ability to 
model coherent carrier transport experiments in low-dimensional semiconductor systems is 
essential for designing possible future devices for coherent electronics or quantum computing. 
In the experiments of Refs. Ill and 13, two paths are electrostatically defined in a high-mobility 
AlGaAs 2DEG, within a multi-terminal setup that allows to overcome the phase-rigidity 
constraint 3 of a two-terminal one. Two narrowings along one of the paths define a QD which 
is operated in the Coulomb blockade regime (a different set of experiments, performed in the 
Kondo regime, presents another peculiar phase behavior-^). The transmission phase across 
the QD is measured by an electron interferometry technique in which electrons are emitted 
at a given energy from a quantum point contact at one end of the two-path system: When 
the energy corresponds to a quasi-bound level (QBL) of the QD, a transmission resonance 
occurs. The depth of the QD confining potential Vd is tuned by charging a nearby "plunger" 
gate and the transmittance, together with the corresponding phase, is obtained as a function 
ofV d . 

The process of electron scattering through the QD has been modeled by means of a num- 
ber of different approaches, ranging from multi-particle few-sites^ to lattice^ and Hubbard 8 
model Hamiltonians 9 -'--. Still, none of the proposed approaches has been able to fully repro- 
duce the main feature of the measured transmission phase 9, namely, the recurring behavior 
found in the many-particle regime of the QD, where 9 smoothly changes by it on each trans- 
mission peak of the iV-electron system, and then abruptly drops to the initial value in each 
valley between the N and N + l resonances, this leading to in-phase transmission resonances. 
This is called the universal behavior since it does not depend on the charge status of the 
QD. While the change of the phase at each resonance is well described by the Breit-Wigner 
model, the nature of the phase drops remains substantially unexplained. 

Recently, an enhanced version of the electron interferometer system-, allowing for the 
precise control of the number of electrons inside the QD down to zero, has been used to 
measure the coherent transmission amplitude for small N. The results show that when only 
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a few electrons (N < 6) are bound into the QD, the universal behavior of the phase is lost, 
and the phase drop occurs only for certain values of N. Furthermore, it was confirmed that 
the measured phase evolution is indeed related to the iV-electron dot and not to the larger 
two-path device. 

The aim of the present paper is to show that the universal behavior of the phase (large 
N) is consistent with an electrostatic approximation, where the electron-electron interaction 
between the scattered carrier and the bound ones is included as a mean Coulomb field. 
This is done in Sec. [TXJ, where the transmission probability and phase are computed for 
a 2D potential representing the QD (attached to source and drain leads) plus a "large" 
number of bound electrons. Furthermore, in Sec. II III we show that an exact few-particle 
calculation performed on an effective ID model of the system leads to the appearance of both 
Breit-Wigner and Fano resonances^, with continuous and discontinuous phase evolution, 
respectively, consistent with the experimental findings in the small N regime. Finally, in 
Sec. HVl we draw our conclusions. 



II. MEAN-FIELD APPROACH: RECURRING PHASE DROPS 

Let us resume the expected phase evolution for a single electron crossing an empty QD. 
We do so for a specific 2D potential V s [Figs. [H(a) andQJb)] which mimics the one generated 
by the surface metallic gates in the 2DEG of the devices of Refs. Ill and fl Along the 
propagation direction [Fig. QJb)], we take two smoothed barriers (with a maximum height 
of 10 meV and a maximum width of 10 nm) that connect a 60 nm flat negative region which 
mimics the QD potential which is tuned in the simulations; along the transverse direction, 
we consider a harmonic confinement with fku =30 meV. In the setup used in Refs. Hand [J no 
bias is applied between the QD source and drain leads since the coherent electron traversing 
the dot is emitted by a quantum point contact at one end of the two-path system (not 
included in our simulations). Accordingly, we keep the Fermi energy of the two leads at zero 
potential and fix the energy of the incoming electron. Material parameters for GaAs have 
been used. The open-boundary single-particle 2D Schrodinger equation has been solved by 
using the quantum transmitting boundary method 12 in a finite-difference scheme. Figure [Tfc) 
shows the transmission probability and phase as a function of the QD potential. As the QD 
potential is varied, the incoming electron comes into resonance with higher single-particle 
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position (nm) dot potential (meV) 

FIG. 1: (color online). The adopted 2D potential profile (a) consisting of a harmonic potential 
with level spacing of 30 meV in the transverse direction and a double barrier along the propagation 
direction (b). Since no bias is applied the Fermi levels in the source and the drain coincide, and are 
taken to be zero. Single-particle transmission probability and phase (c) are shown, as a function 
of the QD potential. The energy of the incoming carriers is 1 meV. 

QBLs. At each resonance peak, the phase increases by ir in agreement with the Breit-Wigner 
model, while it is substantially constant in the low-transmission valleys. 

We show next that when the mean Coulomb field of electrons that populate the QD is 
taken into account, the behavior of the transmission phase shows the observed drops. In 
our model, the maximum of a transmission peak corresponds to the alignment of the energy 
of the scattered electron with the energy of a QBL of the mean-field potential. When the 
energy bottom of the QD is further lowered and the alignment is lost, the transmission 
probability decreases until the QBL becomes a genuine bound state, i.e., its energy falls 
below the Fermi energy, and it is occupied by an additional electron. The new mean-field 
potential has the QBL of the previous resonance shifted by the addition energy and, after 
a further lowering of the QD potential, it produces another resonance. This phenomenon, 
that is essentially a Coulomb blockade effect, is repeated each time a carrier is added to the 
QD. As the mean fields produced by iV or iV + 1 electrons are very similar in the large N 
regime, the QBL that generates the resonances and the corresponding transmission phase 
is always the same at each peak, with an abrupt drop each time a new electron occupies a 
bound state of the QD. 

We now apply our model to a QD with the structure potential V s of Figs. [T](a) and 
Hfb). In order to estimate the QD electrostatic potential we first solve the closed-boundary 
Schrodinger equation then add the field generated by an electron in the ground state ipi, 
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namely 

with r = {x — x') 2 + (y — y') 2 + {d/2) 2 and where d = 1 nm represents the thickness of 
the 2DEG and A^> = 30 nm is the Debye length-^. The Fermi levels of the source and drain 
leads are fixed, i.e., we neglect the effect of the charge inside the QD on the leads. We 
compute the ground state of the new potential V s + V\ and we repeat the whole procedure 
until we reach a number N of bound particles for which the potential V s + V\ + ■ ■ ■ + Vjy has 
an unbound (positive energy) ground state. Then we compute the 2D scattering state for 
an incoming electron with the boundary conditions already described for the single-particle 
calculation. For simplicity the bound states are calculated in a finite domain by solving 
the closed-boundary Schrodinger equation. This leads to a shift in the energy of the bound 
states that has no effect on the qualitative results of the present work, i.e., the phase drops 
between the transmission resonances. 

We show two sets of calculations in Fig. [2J In the top panel (a), the system parameters 
are chosen as in Fig. [1] in order to obtain a clear resolution of the resonances, although 
they do not correspond to the experiments in Refs. Q and 13. For the chosen parameters, 
the transmission occurs through the fourth excited QBL. All resonances, corresponding to 
different N, are in phase and this trend continues as the potential of the QD deepens, i.e., for 
larger numbers of bound electrons. Note that, although the effect of the charging of the QD 
is essentially classical, the transmitted electron must be obviously modeled in a quantum 
approach in order to obtain the transmission phase. 

In Fig. Mjo), we consider a structure with parameters closer to those of the experimental 
conditions in Ref. y; in particular, the confinement potential is much weaker and the energy 
of the incoming electron smaller than in Fig. [2](a) (see caption), leading to less defined 
resonances. The transmission phase evolution is similar to the previous case, in spite of the 
fact that differences between the two calculations are not only quantitative, showing that 
the results obtained are robust against the details of the calculation and of the system. In 
particular, (1) due to the low energy of the incoming carrier, the transmission takes place 
through the ground QBL rather than an excited state; (2) the two lowest QBLs are, for 
N > 4, quasi-degenerate.— The latter effect is due to charge accumulation in the center of 
the QD, away from the barriers, inducing a double-well-like profile along the propagation 
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FIG. 2: (color online), (a) Transmission probability and phase for an electron scattered by the 
potential corresponding to the sum of the structure potential described in Fig. Q] and the mean 
field generated by N bound electrons. The numbers indicate the value of N at each transmission 
peak: as the QD potential Vd decreases iV increases, (b) Transmission probability and phase are 



shown for a system simi 
the energy levels of Ref. 



ar to the one of Fig. [T] whose parameters are tuned in order to match 
2j, namely: kinetic energy of the scattered electron ~ 20 /ueV, charging 
energy ~ 1 meV, coupling of the QD with the leads ~ 200 ^eV, and difference between the first two 
single-particle QBLs ~ 500 fieV. The above values are obtained with a 100 nm well and two 50 fj,eV 
barriers 4 nm wide in the longitudinal direction and a harmonic confinement with ftio = 1 meV in 
the transverse direction. 

direction. In this regime, the transmission peaks corresponding to the two lowest QBLs 
merge and, for each N, a single transmission resonance is found that, being originated by 
two quasi-degenerate states, is characterized by a phase change of 2ir. However, since the 
trapping of an additional electron in a localized state takes place just after the transmission 
maximum, the resulting phase evolution spans only a range of tt. A further effect of the 
charge accumulation in the center of the QD is the decrease of the maximum value of the 
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transmission probability on the resonances. The above trend is clear in the left part of 
Fig. [2](b). We note that our simulations are performed at zero temperature and with an 
exact energy of the incoming carriers, this leading to the steep transitions in the transmission 
probability of Fig.^b). Such steepness is not expected in experiments due to the uncertainty 
of the incoming carriers' kinetic energy and the temperature dependence of bound levels' 
occupancy. In the simulations based on the mean-field approach, the universal behavior, i.e. 
the phase drops occurringbetween successive resonances, persists down to iV = 0, in contrast 
with experiments of Ref. |2| where phase drops may or may not occur for iV < 6. It should be 
noted, however, that the phase drops are a necessary consequence of the electrostatic model 
employed for the coupling between the bound and incoming electrons: Such a mean-field 
picture is expected to break down at small N. Indeed, we show in the following section that 
the inclusion of carrier correlation may give rise to an iV-dependent phase evolution. 

To conclude our mean-field analysis we discuss the similarities between the results hitherto 
presented and the ones obtained in Ref. [s, also including electron-electron interaction in a 
mean-field approximation. In the above work, the lead-dot-lead system is modeled with a 
cross-bar geometry and the transmission amplitude is obtained by using the non-equilibrium 
Green function approach and a Hubbard Hamiltonian. The recurring phase drops are found 
at zeros of the transmission 15 and persist when the electron-electron interaction is turned off. 
While the first effect agrees with our simulations, we find no drops in the non-interacting case. 
The difference can be explained by the different models adopted for the dot: a 2D doub 
barrier structure in our case and a ID bar orthogonal to the lead-to-lead direction in Ref. 
This is confirmed by the further agreement between our Fig. []Tc) and a non-interacting 
simulation for a double-barrier ID structure reported in the above work. There, the universal 
behavior of the transmission phase seems induced by the cross-bar configuration, with a 
single site between the two leads, regardless of the Coulomb interaction. 



III. FEW-PARTICLE APPROACH: BREIT-WIGNER AND FANO RESO- 
NANCES 

In order to obtain numerically the transmission coefficient of a fully correlated system 
the calculation must be able to solve the few-particle problem exactly in an open domain, 
a difficult task for a general 2D potential. We therefore chose to simulate the dynamics for 
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two electrons in a strictly ID quantum wire with the same profile of Fig. E^b); the lateral 
extension of the wire is, however, taken into account by an effective Coulomb potential 
Vc(x) = e 2 /[47re (x + d)}, where the Coulomb singularity is smoothed by a cutoff d — 1 nm^. 
Then, we solve exactly the few-particle open-boundary Schrodinger equation in the real 
space, using a generalization of the quantum transmitting boundary method mentioned 
above, whose general derivation is detailed elsewhere^ 1 ^. In the following, we describe it 
for a ID spinless system. 

Let us consider a region of length L, with a single-particle potential V(x), constant 
outside that region (leads): V(x) = V(0) if x < and V(x) = V(L) if x > L. Although the 
method is valid for the general case we consider here V(0) = V(L) = for simplicity. Let 
us take (N — 1) interacting identical particles bound by V(x) in its (N — l)-particle ground 
state Xo{ x h ■■■i x N-i)- The m-th excited eigenstates of (N — 1) interacting particles will be 
denoted by Xm- 

Our aim is to find the correlated scattering state of iV-particles i[>(xi, . . . ,xn) that has 
the following form when the n-th (with n < N) particle is localized in the left lead, i.e., 
when x n < 0: 



^(x 1 ,...,x n ,...,x N )\ Xn<0 = (-iy 



Xo(xi, x n ~i,x n+1 , x N )e tk ° Xn + 



+ Yl b rnXni( X li X n-1, X n+1 , X N )e~ ik ™ Xn + (2) 
m=0 
M 

+ b mXm(xi, ...,x n _i,x n+ i, ...,x N )e k ™ x 

m=Mi+l 



where (—1)™ accounts for the wave function antisymmetry and k l m = ^J\2m e T m \ represents 
the wave vector of the traveling particle, with mass m e , whose kinetic energy T m = E — E m 
is obtained from the total energy E and the energies E m of the states Xm] i n turn E can be 
obtained from E = Tq + Eq since the incoming-particle energy To is the Fermi energy in the 
left lead. 

The first term inside the square brackets on the r.h.s. of Eq. (j2J) accounts for the n-th 
particle incoming as a plane wave with To energy from the left lead, while the other (N — 1) 
particles are in the ground state of V. The second term represents the linear combination 
of all the energy-allowed possibilities with the n-th particle reflected back as a plane wave 
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in the left lead with energy T m and the dot in the Xm state. The third term is analogous to 
the latter but accounts for the case E m > E, representing the n-th particle as an evanescent 
wave in the left lead. The number of bound states whose energy is lower than the total 
energy, E m < E, is Mi + 1. When particle n is in the right lead, the wave function has a 
form similar to Eq. ([2]), without the incoming term since we are considering only electrons 
traversing the dot from the left lead: 

|- M r 

ip(x 1 ,...,x n ,...,x N )\ Xn>L = (-1)" ^2b^ l Xm(xu...,x n _ u x n+1 ,...,x N )e lk ' mXn 

-m=0 

M 1 

+ b mXm{xi, ...,x n _i,x n+ i, ...,x N )e~ k ™ Xn . (3) 

m=M r +l - 



Since the number of interacting particles is iV and the problem is ID, our computational 
domain consists of an iV-dimensional hypercube that we discretize with a real-space square 
mesh. On the internal points, the wave function ip{x\, . . . ,Xn) satisfies the usual iV-body 
Schrodinger equation 



ft 2 A <r 2 v v " 



2m e ^— ' dxi, 

e n =l n n=l n=l n'=l 



il){x u . . . , X N ) = E ip(x-i_, ...,X N ) (4) 



where V s (x) is the ID single-particle potential energy of the structure at position x and 
Vc(d) is the effective ID Coulomb energy of two electrons at a distance d. E is the total 
energy, defined previously. 

The wave function ip has to match Eqs. ([2]) and ([3]) on the N "left" boundaries and the 
iV "right" boundaries of the domain, respectively. Equations (J2J), ([3]) and (jl]) are then solved 
together as a coupled system of 2N + 1 equations, using a finite-difference discretization for 
the derivatives. In this way, the reflection amplitudes b < and the transmission amplitudes 
6 > (i.e. the unknowns, together with ip) are obtained numerically— Note that the result- 
ing wave function is antisymmetric for any two-particle exchange since we have imposed 
antisymmetric boundary conditions. 

Figure [3]^a) shows the transmission amplitude for a two-particle correlated triplet state 
(bound and traveling electrons with the same spin) as a function of the QD potential. 
The energy of the incoming electron is T = 20 /xeV, as in the mean-field simulations of 
Sec. [Til The first five transmission resonances are clearly visible, three of which have a 
Lorentzian shape with a width of about 200 /ieV and a Breit-Wigner-type phase evolution, 
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FIG. 3: (color online), (a) Correlated transmission probability (bottom) and phase (top) for an 
electron scattered by the potential described in Fig. [2] caption when a second electron is bound in 
the QD. Three Breit-Wigner and two narrow Fano resonances are present, (b) and(c) Details of the 
transmission spectrum, with transmission probability in logarithmic scale, showing the asymmetric 
Fano resonances and the corresponding phase jump of tt. 

similar to the ones already seen in Fig. [TJ The two remaining resonances, shown in detail 
in Figs. [3]^b) and[3](c) (note the logarithmic scale for the transmission amplitude), are very 
narrow (few fieV) and present a typical asymmetric Fano line shaped. They are a signature 
of electron-electron correlation and, from their small width, we deduce that the effect of the 
Coulomb potential is very limited in our model. Nevertheless, the distinctive behavior of 
the transmission phase is clearly visible in the upper plots of Figs. [3](b) and[3](c). An abrupt 
phase jump of 7r takes place near the resonance, where the transmission probability vanishes. 
This shows that the origin of the phase jumps detected in the few-particle experiments may 
reside in correlation-induced Fano resonances. A similar behavior, with the presence of both 
Breit-Wigner and Fano resonances, showing continuous and discontinuous phase evolutions 
respectively, is found in the simulation of the correlated three-electron scattering state (not 
shown here). In the latter case, the ratio of Fano resonances is larger and keeps increasing 
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with the number of particles. 



IV. CONCLUSIONS 

In summary, we showed that in the few-particle correlated regime, both Breit-Wigner 
and Fano resonances are found, while in the mean-field regime, a single type of resonance 
is present, which is repeated for each number of bound electrons, this leading to a phase 
jump of 7r each time a new electron enters the QD. While these results are consistent with 
experiments, the microscopic nature of the recurring phase drops in the latter regime remains 
unclear. In fact, they are well reproduced by the first of our approaches in which no quantum 
correlation is present between the bound electrons and the scattered one. On the other hand 
different models 1 - that take into account the fully-correlated dynamics of the carrier such 
as our few-particle calculations cannot be applied to the many-particle regime. The above 
considerations suggest that in the latter regime, the coherent component of the transmitted 
electron wave function (i.e., the component that does not get entangled with the QD and 
whose phase is detected by the interferometer) behaves as if the QD, together with the bound 
electrons, was a static electric field, being unable to discriminate between two different values 
of N. On the other hand, when N is small, the transmission phase is able to provide a 
partial information on the number of electrons confined in the QD through the character 
of the transmission resonance and the possible phase lapse. The transition regime between 
the many- and few-particle conditions" needs further analysis since it can clarify the 
connections between the two opposite approaches used in the present work. 

We finally note that the Fano resonances found in the ID two-particle scattering states are 
a genuine effect of carrier-carrier correlation, in accordance to the original concept developed 
in Ref. |2C|. A more general definition is often adopted, ascribing the Fano line shape of the 
transmittance to the interference of two alternative real-space pathways. In fact, Fano 
resonances have been obtained previously by means of 2D 23 , multi-channel 1 ^, and two- 
path' 1 ^ single-particle calculations. In those cases, however, the ratio between the number 
of Fano and Breit-Wigner resonances is not expected to vary when varying the confinement 
energy of the QD, while the number of correlation- induced Fano resonances becomes shortly 
dominant when moving from a low- A?" to a high-iV condition in the framework of a full 
few-particle modeling. 
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Upon completion of this work we learned about a recent work by Karrasch et al.— where 
the ii lapses are also ascribed to (anti)resonances of Fano type. 
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